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Electron-phonon coupling in semiconductor quantum dots plays a significant role in determining 
the optical properties of excited excitons, especially the spectral nature of emitted photons. This 
paper presents a comprehensive theory and analysis of emission spectra from artificial atoms or 
quantum dots coupled to structured photon reservoirs and acoustic phonons, when excited with 
incoherent pump fields. As specific examples of structured reservoirs, we chose a Lorentzian cavity 
and a coupled cavity waveguide, which are of current experimental interest. For the case of optical 
cavities, we directly compare and contrast the spectra from three distinct theoretical approaches to 
treat electron-phonon coupling, including a Markovian polaron master equation, a non-Markovian 
phonon correlation expansion technique and a semiclassical linear susceptibility approach, and we 
point out the limitations of these models. For the cavity-QED polaron master equation, we give 
closed form analytical solutions to the phonon-assisted scattering rates in the weak excitation ap¬ 
proximation, fully accounting for temperature, cavity-exciton detuning and cavity dot coupling. We 
show explicitly why the semiclassical linear susceptibility approach fails to correctly account for 
phonon-mediated cavity feeding. For weakly coupled cavities, we calculate the optical spectra using 
a more general reservoir approach and explain its differences from the above approaches in the low 
Q limit of a Lorentzian cavity. We subsequently use this general reservoir to calculate the emission 
spectra from quantum dots coupled to slow-light photonic crystal waveguides, which demonstrate a 
number of striking photon-phonon coupling effects. Our quantum theory can be applied to a wide 
range of photonic structures including photonic molecules and coupled-cavity waveguide systems. 

PACS numbers: 42.50.-p, 42.50.Ct, 42.50.Nn, 78.67.Hc 


I. INTRODUCTION 

Artificial atoms or quantum dots (QD) are excellent 
candidates for solid-state quantum bits (qubits) and show 
promise for enabling scalable quantum information pro¬ 
cessing at optical frequenciesi’-. Embedding QDs in pho¬ 
ton cavity structures can manipulate the radiative de¬ 
cay rate by tailoring the local photon density of states 
(LDOS). For example, one can enhance the spontaneous 
emission (SE) rate through the Purcell effect^, or sup¬ 
press SE leading to long lifetimes of around tens of ns^, 
which is promising for low error rate quantum logic oper- 
ations^. However being part of a solid-state lattice struc¬ 
ture, QD excitons are intrinsically coupled to the under¬ 
lying phonon bath-, which significantly reduces their co¬ 
herence time on short time scales. Phonon dressing of QD 
emission manifests itself in a number of experimental ob¬ 
servations, such as phonon-assisted inversion-!”— , damp¬ 
ing and frequency shifts of driven Rabi oscillations!^”— 
and excitation induced dephasing of Mollow triplet side- 
bands!!’-!^, which distinguishes QDs from simple two- 
level atoms!2. Over the past decade, several theories 
have been developed to address this issue. They range 
from the independent Boson model (IBM)^’— , correla¬ 
tion expansion techniques!!, perturbative master equa¬ 
tions (MEs)— , polaron MEs22 — , variational MEs2! and 
path integral calculations!!. 

Coupling QDs to structured photonic reservoirs like 
photonic crystals can be used for building a scalable 
quantum light-matter interface! and allows quantum 
control of light-matter interactions!. One common de- 



FIG. 1: (Color online). Two selected schematics showing a QD 
in a structured photonic reservoir. Photonic microcavity (a) and 
a coupled-cavity waveguide (b) using a photonic crystal platform, 
containing an artificial atom or a coupled semiconductor QD. 


scription for the QD-medium light interactions is to treat 
the QD like a simple atom with discrete energy levels. 
However, the phonon interactions may strongly influ¬ 
ence the light-matter coupling and in turn modify the 
emission properties of the QD in a non-trivial man- 
ner-!!”-!!. For example, the well known Fermi’s golden 
rule breaks down in weakly coupled photonic reservoirs 
which have photon correlations times that are compa¬ 
rable the the phonon bath and this results in enhanced 
or even suppressed SE!!. Experimentally, phonon inter¬ 
actions are clearly visible in the photon emission spec¬ 
tra, especially when coupled to a high Q cavity, where 
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Q is the quality factor. When electron-phonon coupling 
is included in the models, polaronic approaches'^ and 
linear susceptibility theories^’— (where the IBM model 
is included in the QD polarizability) predict asymmet¬ 
ric vacuum Rabi splitting—’—’— and strong off-resonant 
cavity feeding—’—’—, which have been observed experi¬ 
mentally—-—. 

In this paper, we present a comprehensive analysis of 
emission spectra from QDs coupled to structured pho¬ 
tonic reservoirs, in the presence of electron-phonons scat¬ 
tering. The development of a general light-matter inter¬ 
action theory for a QD exciton coupled to an arbitrary 
photon and phonon bath is extremely challenging and 
several approximate treatments have been employed to 
treat several regimes of interest. For example, the cavity- 
QED (cQED) polaron ME approach, which treats the 
cavity mode operator at the level of a system operator, 
is suitable for modelling high-Q cavities^. If a weak- 
to-intermediate coupling regime with the photon bath 
is assumed, the polaronic ME theory can be extended 
to handle arbitrary structured photon reservoirs using a 
photon bath approximation—. A semiclassical approach 
has also been used to model the spectrum in the linear 
excitation regime22’.25^ though is restricted in its range of 
excitation conditions and modelling true quantum light- 
matter interactions; we show explicitly the failure of such 
an approach for modelling cavity feeding over large QD- 
cavity detunings. Without cavity and photon bath cou¬ 
pling, a common perturbative approach to describing 
electron-phonon coupling has been through the correla¬ 
tion expansion technique, which has been used to ex¬ 
plain phonon-induced damping of Rabi oscillations where 
electron-phonon coupling is treated to second-order—’—; 
we extend such an approach to study QDs in a simple 
cavity and make a direct comparison with the polaron 
ME approaches and the linear susceptibility model. 

Each of these approximate theories has its own mer¬ 
its and limitations, though there has been little investi¬ 
gation into how each of these approaches compare with 
each other, and, to the best of our knowledge, there has 
been no quantum approach for treating the emitted spec¬ 
trum for excited QDs in the presence of phonons and 
general photon reservoirs. Eor the photon medium, we 
consider examples of a simple photonic crystal cavity-1 
(Eig. l.(a)) and a coupled resonator optical waveguide 
(CROW)^ (Fig- 1(b)) which are relevant in the con¬ 
text of current semiconductor QD experiments. Eor the 
case of a high-Q cavity, we examine the case of phonon- 
dressed vacuum Rabi sphtting-1^’— , when strongly cou¬ 
pled QD-cavity are at resonance. Eor the far off-resonant 
case between the QD exciton and a Lorentzian cavity, 
we consider the scenario of phonon-assisted cavity feed- 
ing-11 ’— Pqj. photonic cavity, the polaron- 
transformed QD-cavity-H interaction is simplified using 
a Born-Markov approximation-11’—. This approximation 
is expected to be valid if the system dynamics is longer 
than the correlation time of the phonon bath. To investi¬ 
gate the possible breakdown of this approach, we utilize 


the non-Markovian correlation expansion approach-H to 
include a cavity, and we use this method to help dis¬ 
tinguish the regimes in which the polaron ME approach 
is valid, since the latter is considerably simpler to use 
and can easily include in other quantum processes such 
as multi-photon effects. We also compare the above two 
approaches with the predictions of a semiclassical but 
non-Markovian linear susceptibility theory-H’-H, where 
the IBM lineshape is incorporated within the frequency- 
dependent polarizability function of the QD exciton. The 
IBM lineshape shows acoustic phonon sidebands that sur¬ 
round the zero phonon line (ZPL). We show the clear 
breakdown of this approach and introduce a fully quan¬ 
tum approach to treat the emission spectrum of a QD 
in the presence of any spectral shaped photon reservoir 
coupling-^, which is valid in the Purcell regime, i.e., for 
weak-to-intermediate coupling. In the case of a simple 
leaky cavity, differences from the cQED polaronic ME 
and the reservoir polaron ME approaches arise, when the 
coupled cavity becomes increasingly lossy-^. Lastly, we 
use this general reservoir approach to investigate emis¬ 
sion spectra of QDs from photonic crystal coupled cavity 
waveguides22. To the best of our knowledge, this is the 
first theory to study such spectral features from waveg¬ 
uides, without making a mean-field approximation for 
phonon coupling effects^. In the mean-field limit, it has 
been common to assume that phonons reduce SE in a 
structured reservoir by a frequency independent constant 
where (B) is the the thermal average of phonon 
bath displacement operators B^— (also see Sec. II A) 
which depends on the bath temperature. The polaronic 
reservoir approach-22, however, shows that SE can also be 
enhanced by phonons and depends on the local density 
of states (LDOS) of the photon bath reservoir; also note 
the polaronic approaches include effects that are nonper- 
turbative in phonon electron-coupling and can be used 
to model low to high temperatures of the phonon bath. 

Our paper is organized as follows. In Sec. II, we intro¬ 
duce the key theoretical models. Sections IIA and IIB 
briefly review the polaron reservoir ME22 and the po¬ 
laron cQED MEi22, respectively. Section IIC presents the 
correlation expansion approachi^ with the coupling to a 
leaky cavity mode. In Sec. IIP, we derive the incoherent 
emission spectrum when the QD is excited by a weak in¬ 
coherent pump field and show how this is obtained from 
the two-time correlation function; and we also present 
the emitted spectrum from an excited exciton using lin¬ 
ear susceptibility theory (Sec. IIP4b In Sec. HE we dis¬ 
cuss the phonon parameters used for calculations, and in 
Sec. HE we discuss the numerical complexity of the dif¬ 
ferent approaches. In Section III, we present and discuss 
various numerical results. We begin with a discussion 
of emission spectra from a bare (i.e., with no photonic 
bath coupling) QD (Sec. Ill Ai and thereafter present re¬ 
sults on QDs coupled to structured photonic reservoirs. 
Section III B investigates the scenario of phonon-dressed, 
resonant vacuum Rabi splitting (i.e., in the strong cou¬ 
pling regime) and reveals the limitations imposed by the 
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Born-Markov approximation on existing polaron cQED 
approaches. Section III C treats the case of off-resonant 
cavity feeding under weak coupling conditions. The ef¬ 
fects of accounting for intercorrelated photon and phonon 
bath dynamics is discussed for low Q cavities using the 
photonic reservoir theory— (Sec. Ill C lb The differ¬ 
ent predictions of the linear susceptibility theory in the 
context of off-resonant cavity feeding is also explained 
(Sec. Ill C 2b Finally, in Sec. Ill D. we study in detail the 
emission spectrum from QDs weakly coupled to at pho¬ 
tonic crystal waveguides; we use a model LDOS suitable 
for a slow-light CROW. Section IV summarizes the main 
results and discusses future directions. In Appendix A 
we show the equivalence between the spectra obtained 
from a weak incoherent pump and an inverted exciton 
as the initial condition, since these two approaches are 
used with the MEs and linear susceptibility model, re¬ 
spectively. 


II. THE HAMILTONIAN WITH PHOTON 
BATH AND PHONON BATH COUPLING TO A 
QUANTUM DOT EXCITON 

We consider a single neutral QD exciton (strong con¬ 
finement limit) that is modelled as a two-level system 
coupled to a photon bath and a phonon bath described 
by lowering operators f and 6^, respectively (see Fig. 2.). 
Such a two-level approximation is valid for small epitax¬ 
ial QDs (e.g., self assembled InGaAs/GaAs QDs), whose 
parameters are used for the current study. Although 
real QDs have many exciton levels over a broad band 
of frequencies, the coupling to one exciton has success¬ 
fully explained a number of experiments when probing 
single exciton dynamics, e.g.. Refs. 15 in addition, 
one can extend such an approach to add in other exciton 
levels^. The deformation potential coupling with longi¬ 
tudinal acoustic (LA) phonons play the strongest role in 
such QDs. The resultant phonon side-bands span a fre¬ 
quency range of ^ ± 5 meV around the s-shell transition 
which is much smaller than the energetic separation be¬ 
tween 5-shell and p-shell transitions (~ 25 meV)i2, thus 
validating a two-level approximation. The lowering oper¬ 
ator a~ describes a transition between the QD states |e) 
and 1^), separated in energy ujx. In a frame rotating at 
the frequency of the QD, the coupled system is described 
by the following Hamiltonian—: 
poc 

dr / dco’f^(r,cc;,t)f(r,CL;,t) + 

q 

duo d • E^(r(^, cj, t) H- H.c. 

+ a+a~'^hXq{bl + bq), ( 1 ) 

q 

where a dipole and rotating wave approximation is made 
to describe coupling between the QD and the photon 
bath. In (A), a QD of dipole moment d is assumed to 
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H = h 


/ 



FIG. 2: (Color online). Energy level diagram of a neutral QD 
(electron-hole pair) interacting with a phonon bath and a photon 
bath. The operator fj^ (bq) creates a photon (phonon). 


be located at the spatial position and the exciton- 
phonon coupling strength Xq is assumed to be real—. 
The positive frequency component of the electric field 
operator E+(r(^, cj, t) can be expressed in terms of the 
electric field Green function G(r,r';co’) as E+(r,co’,t) = 

i f dr'G(r, r'; co’)y^^^^ei(r', cJ)f(r', cu, t)— ; in a medium 

described by the dielectric constant e(r,CL;) = eR(r,C(;) -h 
iei(r,cc;), G(r, r') is the solution to Maxwell equations at 
r to a point dipole oscillating at r' (without any QD cou¬ 
pling). This expression for the electric field operator is 
quite general and fully satisfies the Kramers Kronig rela¬ 
tions in a general photonic medium. In order to avoid any 
coherent pump induced dressing of QD emission (which 
is interesting but beyond the scope of the present paper), 
the QD is assumed to be weakly excited by an incoher¬ 
ent pump P— , which maps on to a range of typical ex¬ 
periments that are performed to measure the low pump 
emitted spectrum for a QD-photonic structure. In the 
linear regime, one can obtain exactly the same expression 
starting from an excited exciton, but this requires one to 
integrate over two times instead on one; this equivalence 
is shown explicitly in Appendix A. In a ME approach, the 
incoherent pumping term is included as a Lindblad oper¬ 
ator, ^L(cr+), where L(0) = 20p0^ — 0^0p — pO^O and 
P is the pump rate; the ME can also include other Lind¬ 
blad operators to account for additional incoherent pro¬ 
cesses such as background spontaneous decay, ^L(cr“), 
and QD pure dephasing, Neglecting cou¬ 

pling to a structured photonic reservoir and phonons, 
these incoherent processes determine the line-shape of 
the QD emission spectrum, which is represented by a 
simple Lorentzian broadening of the ZPL. 

We stress that the above Hamiltonian is applicable for 
arbitrary photon and phonon baths and is very difficult to 
solve in general. In the following subsections, the Hamil¬ 
tonian is thus manipulated using different approxima¬ 
tions which includes photon coupling, phonon coupling, 
and decoherence processes in a self-consistent and phys¬ 
ically meaningful way. We then use these approaches 
to compute the linear spectrum in different QD photon 
reservoir coupling regimes, and we explore the range of 
validity of the various models. 
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A. Photon-reservoir polaron master equation 
approach 


We begin with the Hamiltonian in (JL) and perform a 
polaron transformation. This unitary transform includes 
phonons to all orders but puts the Hamiltonian in an eas¬ 
ier form to apply system-bath reservoir theory and per¬ 
turbation theory, with polaron shifted interaction terms. 
Essentially we are employing a basis change using a uni¬ 
tary transform in which the bath-shifted exciton becomes 
the polaronic quasiparticle. In this way, we fully recover 
the IBM without photon interactions and can include 
processes that give rise to the ZPL. 


The polaron transformation is given by H' 

where P = , which 

yields a polaron-transformed Hamiltonian, 


dcc;f^(r,cc;, t)f(r,co’, t) -h hujqblbq, 


pOC 

/ dcj d • E^(rc^, cj, t) + H.c. 
Jo 


( 2 ) 


where B± = exp[±Y,q^{bq - &J)] are the coherent 
phonon bath displacement operators22. The polaron shift 
Ap = duj also appears in the polaron trans¬ 

form and below we will assume this factor is absorbed 
in the polaron-shifted frequency of the QD, defined by 
{= ujx — Ap). In the continuum limit, the phonon 
coupling is determined by the phonon spectral function 
• If a weak interaction between the QD and the 
photon reservoir is assumed, a ME for the QD reduced 
density matrix p can be derived by retaining terms up to 
second order in the polaron-shifted interaction Hamilto¬ 
nian, H[ = —[Bpcr^e^^Pj^dujd • E+(r(^,cj, t) + H.c.], 
which is also formally known as the Born approximation. 
The form of the time convolutionless^ (or time-local) 
interaction picture ME is then 


dp{t) ^ 
dt 

(3) 


where H'i{t) = ex.p[iH^t/h]H[ex.p[—iH^t/h] and = 
hf dr duo f^(r, uo, t)f(r, uo, t) P hwqhyoq. The trace 
operators TrR^j^ and TrRp^ denote a trace over the photon 
and phonon reservoirs, which are assumed to be statis¬ 
tically independent with pR = PRphPRpn— • The trace 
over the photon reservoir—’— assumes thermal equilib¬ 
rium and uses the relations, TrR^^ [f(r, cj), f^(r', cj')] = 
[n(w) + l](5(r-r')<5(w-w') and TrRpj,[r(r,w),f(r',w')] = 
h{uo)d(y — y')5{uo — uo') where n{uo) ~ 0 , which is valid for 
optical frequencies. The final form of the polaron reser¬ 
voir ME in the Schrodinger picture^! and Markov limit 


(t ^ oo in the integral of (^)) is 


^ ) - *ALamb [o-+cr , p] + yi(cr ) 

+ :^£(o-+cr") + ^i(cr+), (4) 


where the SE rate of the QD into the structured reser¬ 
voir—’— is 


7 = 2/ Re[C'p„(r) Jph(r)]dT, (5) 

Jo 

and the QD Lamb shift is ALamb = 
Im[Cpn(T) Jph(r)](ir. Note that the QD back¬ 
ground decay, pure dephasing and incoherent pumping 
is also included in the ME. The resulting SE rate and 
the Lamb shift display an interplay between photon and 
phonon bath dynamics, where Jph(^) and Cpn(T) are 
the photon and the phonon bath correlation functions, 
respectively. The photon bath correlation function is 
defined as Jph(^) = j^duoJ^]^{uo)e^^^^~^^'^^ where the 
photon reservoir spectral function is 


Tph(c^) 


d • Im[G(r^,r^;cj)] • d 
irheo 


The phonon correlation function 

Cpnir) 


(6) 

(7) 


where the IBM phase function is defined through 



which includes a sum over phonon emission and aborp- 
tion processes. Note that the SE rate 7 in principle in¬ 
cludes contributions from photonic LDOS values at fre¬ 
quencies different from the ZPL frequency of the QD 
(where uo = uo'^). Such non-local frequency effects for 
the SE is caused by a breakdown of Eermi’s golden rule, 
which depends on the correlation times of both photon 
and phonon reservoirs22. We note that (A) is broadly ap¬ 
plicable irrespective of the specific structure of the pho¬ 
ton reservoir. This approach is restricted to weak-to- 
intermediate coupling between the QD and the photon 
reservoir, and thus it cannot treat the strong coupling 
regime which would manifest in effects such as vacuum 
Rabi spitting; however, this is generally the case for a 
photon “reservoir.” To treat the strong coupling regime 
of quantum optics, e.g., for a high-Q cavity, one includes 
a cavity photon operator at the level of a system operator 
which requires a different approach based on the polaron 
cQED ME that we describe below; while being able to 
describe the strong coupling regime, this approach is re¬ 
stricted (in the presence of phonon coupling) to model 
simple cavity structures and may also breakdown for low 
Q cavities22. 
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B. Cavity-QED polaron master equation theory 


To derive the polaron cQED ME, we replace the pho¬ 
ton bath term in the Hamiltonian given by (JL) with a 
single cavity mode, so that 

H = h/S^cxCi^CL -h hujqb^qbq + ^(cT^a -1- aV“) 

Q 

+ + H.c.) + cr+(7“ ^ + bq), (9) 

m q 

where Ac^c = ooc — ^x is the cavity-QD detuning, g is the 
QD-cavity coupling for a single cavity mode described by 
lowering operator a, and ftrn represents the coupling to 
the photon environment that causes decay of the cavity 
mode. The time-dependence of a is kept implicit in (£). 
The QD-cavity coupling g can be expressed in terms of 
the dipole moment d = dh where g = 77 ( 11 , ^ 5 

where Vqh is effective mode volume of a dielectric cavity 
with dielectric constant e, and 77 accounts for any devi¬ 
ation from the field antinode position and misalignment 
in polarization coupling, and for optimal coupling is sim¬ 
ply 1 . A dipole and a rotating wave approximation is 
used in describing the QD-cavity interaction. As before, 
we perform a polaron transform on H and the polaron 
transformed form— is now given by 

H' = h/S.cx'Ci^CL + hujqb^qbq 

q 

+ g'CiCi^cr )Xg(gXuCui ( 10 ) 

where g' = {B)g^ Xg = g[(j^a + aV“], Xu = ig[(T^a — 
aV“] and the phonon fluctuation operators (^g = -|- 

H_-2 {B)) and Cn = where {B) = (H+) = 

(H_) = The polaron shift Ap is once again ab¬ 

sorbed into the QD frequency and Xcx' = uJc — ^x- 
We can then derive a time-convolutionless ME for the 
reduced density matrix p of the QD-cavity system; fol¬ 
lowing Ref. we use a second-order Born approxi¬ 
mation with the polaron transformed interaction H[ = 
Xg(^g -h XuCu' The master equation in the interaction 
picture is then 

m = dt'TrnmmminmPB]]}, (n) 

where Hi{t) = with = 

HAcx'Cl^cl + g'[a^a + aV“] and hujqb^gbq. The 

operator Trp denotes trace over the phonon bath (pp) 
and we assume the full density operator to be separable at 
all times ppe— • Performing the trace and transforming 
back into the Schrodinger picture, we obtain the final 


from of the time-local Markov polaron cQED ME25, 

+ ) + ^L{ct a ) + -L{ct ) 

'lo m=g,u 

+ H.c.) + |i:(a), ( 12 ) 

where k is the cavity decay rate, and Gg = 
{B)‘^ (cosh((/)(t)) — 1) and Gu = (H)^ sinh(0(t)) are the 
polaron Green functions^S; as before, we have extended 
the upper limit of the integral in (11) to t ^ 00 to ob¬ 
tain a Markov form. It should be noted that the inco¬ 
herent processes are also included in the ME using the 
respective Lindblad terms. Eor the Born-Markov approx¬ 
imation to be valid, in the polaron frame, the system 
dynamics should be much slower compared to the rate 
of relaxation of the phonon bath. Although the bath 
relaxation time is only a few ps for typical QDs, this ap¬ 
proximation may be restrictive in certain regimes, e.g., 
when dealing with vacuum Rabi oscillations at large g 
(Sec. IIIB) and off-resonant cavity feeding at large de¬ 
tunings Aca; (Sec. Ill Cb In order to investigate this possi¬ 
ble limitation, we extend a previously proposed phonon 
correlation expansion approach in Sec. IIC by adding 
a cavity which is not limited by the Born-Markov ap¬ 
proximation, though we only include phonon coupling 
to second-order. Although the correlation expansion in¬ 
cludes phonon effects perturbatively, the approach shows 
good agreement with the full IBM model for linearly ex¬ 
cited bare QDs (i.e., with no photon bath) at low tem- 
peraturesi^. Photon coupling effects with the correla¬ 
tion expansion approach have been highlighted in Ref. 47 . 
though the exciton-photon correlation was truncated us¬ 
ing a Markov approximation, which is strictly valid only 
in weak coupling condition or for unstructured photon 
reservoirs. We highlight two useful limits of the cQED 
polaron ME: when phonons interactions are turned off, 
one fully recover the Jaynes-Cummings model; and when 
the cavity is turned off, one recovers the IBM—, with the 
addition of important ZPL processes. 


Before concluding this section it should be noted that 
(12.) can be further simplified to a simpler effective master 
equation—, if a weak excitation approximation (WEA) is 
made, which is exact for the linear spectrum. The cQED 
polaron ME now takes the analytical form. 


I T ^ ^ I T ^ ^ I ^ 



I/(aV ) 


-h {Mi[(aV“ -|- cr^a), (2cr^cr“a^a + 
+ {M2[(aV“ — cr^a), (2cr^cr“a^a -1- 


— a^a 

— a^a 


)/>]+H.c.} 

)p]+H.c.}, 

(13) 
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where Hf = a+aa^a-+hA‘^^^a'fa-a+a. If 

w e denote the QD-cavity system Rabi frequency as = 

+ 4g'‘^ , then the phonon-mediated cavity/exciton 
scattering rates are given by 

pa^cr“/cr'*”a _ 


■ poo ^ /2 

J ~ cos(lIr))(e“‘^^^^ - 1 ) 


2g'^Re 

+(^ (1 - cos(nr)) + cos(nT)) - 1) 


± 2g'Rm 


/ dT^sin(flT)(e‘^(^) - 1) 

Jo 


(14) 


and the phonon-mediated Lamb shifts are given by 


^a'cr j a 


poo ^ /2 

J ~ cos{^T)){e~'^^'"'> - 1 ) 


^ Im 

+(^ 2 "(l - cos(rir)) + cos(llT))(e'^^'^^ - 1) 


.12 


Tg'^Re 


/ dr ^ sin(Or)(e'^^^^ — 1 ) 

Jo 


(15) 


We also have a cross-dephasing term— 

7cd = 2^'^Re J <^^(-^(1 - cos(flr)) 

-hcos(flr))(e“'^^'^^ ~ 1) + “ cos(flr))(e'^^'^^ — 1) 


/ dr ^ sin(flr)(e~'^^^^ — 1) , 

Jo 


— 2ig'^Re 
and the M terms are defined through 

poo 

Ml = - 25'2 / 

Jo 


(16) 


dr ^ (cos(flr) — l)(cosh0 — 1), 

0 ^ 

(17) 

P^ q' 

M 2 =—2ig'^ / dr^ sin(flr) sinh 0. (18) 

Jo 


A physical understanding of the Mi ^2 scattering terms 
can be obtained by deriving the Bloch equations for a 
and a ~, using the WE A. For example, 


d (a) 
dt 

d{^~) 

dt 


+ 

= -gc(J - i{/^cx' + A^ Ja - 


= —gxd — iA"^^ (J 



(19) 


where = n-\- and E|^ = 7o + 7d + ^ + are 

the effective dephasing rates and gc = ig' — Mi — M 2 and 
dx = ig'+Mi—M 2 are the complex couplings of the cavity 
and QD, respectively, in presence of phonons. Thus the 
processes denoted by the scattering terms Mi^ 2 , result in 
a complex coupling between QD and cavity, in the weak 
excitation approximation. At resonance {Acx' = 0) since 
Ml = 0, so the complex QD-cavity coupling is given by 


9c/x = ig' - M 2 . 

Notably the above ME has solved the incoherent scat¬ 
tering terms exactly and one could use such an approach, 
e.g., to investigate the strongly coupled spectra as a func¬ 
tion of temperature. The Rabi oscillations appearing in 
the integrals ensure that the phonon bath is correctly 
coupled to the dressed resonances of the system. 


In the weak coupling limit, specifically when Acx' ^ 


g'^ the Rabi frequency D 


Acx'-) and the inco¬ 


herent cavity and exciton scattering rates are given 


by, r“ 


/ M a _ 


= 2^'^Re[/Q^ — 1)] and 


phonon modified Lamb shifts are given by, Aq — 

g'‘^lm[J^ — 1)]. The phonon mediated 


scattering rates and Lamb shifts are the same as those 
derived in an earlier work—, where an effective Lindblad 
form for the cavity-QED polaron ME was introduced. 


A connection between the polaron reservoir ME ap¬ 
proach (Sec. II A) and the polaron cQED theory can be 
made by expressing the SE rate (A) in terms of these 
Lindblad decay rates^. 


7 = 7P = rg'--+ 2 / 


^ K+rg _Efl-!—) 


^Ix' + (■ 


■^at CT 


( 20 ) 


which, however, is only valid when the spectral width of 
the cavity is much smaller than the width of the phonon 
bath function (^ 5 meV). This condition is satisfied by 
high Q cavities with /^ < 0.1 meV. Moreover, when k ^ 
r7“ - , then 

* = ( 21 ) 

which can be interpreted as a cavity-feeding term plus a 
phonon-modified (via g ^ g') cavity-induced SE rate. 


C. Correlation expansion approach 


To derive a correlation expansion approach-H’iS, we 
first derive the Heisenberg equations (O = — ■|[0,iE]) 
for a and (t~ from (£), 

d =—iAcxd — igo'~ ^ ( 22 ) 

a-=-iga-i'^Xq{a~bl+a~bq), (23) 

q 

where we again make a WE A, i.e., az = —1, to truncate 
the higher-order photon correlations. Again, this is not 
restrictive in the current situation, since we excite the QD 
weakly to obtain the linear spectrum. The optical Bloch 
equations for (a) and (cr“) are obtained by performing 
ensemble average on the Heisenberg equations of motion. 
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(O) = --i {[0,H]) + Tv[OLp], giving 

(a) = - (^iAcx + I) (a) - ig {a~), 

/ -\ • / \ f P + lO + ld\ / -\ 

(cr ) = -ig (o) - ( - - - 1 {a ) 

-i X] ( 24 ) 

q 

where the Lindblad operators for background radiative 
decay, pure dephasing and incoherent pumping have 
been included. The time-evolution of (cr“) depends on 
phonon-assisted correlations {a~bq), whose equations are 
given by 

{a-bg) = -ig {abg) - (^iug + -^ + 7o +7d ^ {cr~bg) 

2 ^ ^ Xrn bqbfYi') ^ ^ ip 5 

m m 

= -^9 (®^g)c “ (^^^9 + (o'“^g)c 

- i\q{nq + 1 )CT~ 

~ ^ ^ (<^ bqbfn^^ — i ^ ^ Xm (<^ ’ 

m m 

(25) 

hujq 

where riq = _ i)-i ig the average phonon oc¬ 

cupation number for mode q given by the thermal 
equilibrium Bose-Einstein distribution at temperature 
T. The simplification in the second line of (25) is 
performed using the correlation expansion technique^S. 
The operator average can be expanded as 

(cr“6,6t) = {^~) ih) {b}n} + {^~K) (^m) + (o'“^m) {K) + 

(cr“) {bqb}^) + where the three operator av¬ 

erage is separated using cluster expansion techniques 
into singlets ((O)), doublet operator averages {{OO')) 

and a correlated part (OO'O ^ . The pure phonon 
correlations such as and derivatives as (6^), 

blbq^ are neglected, as they represent phonon coher¬ 
ences out of equilibrium and the approximation is sim¬ 
ilar to a bath approximation—. Thus = Uq and 

{^D = i^q) = 0 we have {(j-bqbl^) = (cr“) {bqbl^) 
(p~bqbl^) We next observe that {(J~bq) depends 
on higher-order phonon-assisted correlations such as 
{(J^bqbjn)^ and (p~bqhl^^^ and this leads to a hierarchy 
of coupled equations. In the current work, the hierar¬ 
chy is truncated by dropping all correlations involving 
three or more phonon operators, such as {<J~bqbmbn) 
Thus, using this truncation approximation the equations 
for second-order phonon assisted correlations are given 


by 

(cr bqbjYi)^ = ig {cibqbjYi)^ 

(■( I ^l ^ + 7o+7ciA / -77 \ 

( ^-(cjg “h UJrn) 2 j ' 

i (p bq^^ XYnipm T 1) i {p Xqipq + 1). (26) 

The equations of motion for the phonon-assisted cavity 
correlation are as follows: 

P)c = “ {ii.^cx + W,) + 0 {abq)^ - ig {cr~bq)^ , 

(27) 

{Ojbqbrn)^ — (^(^ca: "t“ "t“ ^m) "t“ '^) {^bgbm}^ 

- ig {(J~ bqbm)^ ■ (28) 

The truncation approximation helps us to derive a 
closed set of equations in the one-time operator averages 
((cr“), (a)), first order ((cr“6g)^, {abq)^) and second-order 
{{(J~bqbm){abqbm)^) phonon correlations, which can be 
solved using direct numerical integration. The emission 
spectrum (see Sec. IIP) requires the two-time operator 
averages or correlation function as (cr+(t)cr“(t -h r)) and 
(^a^{t)a{t + t)). These two-time correlations can be cal¬ 
culated using the one-time coupled equations and the 
quantum regression theorem—. The quantum regres¬ 
sion theorem states that if single time operator averages 
for a set of operators On{t) satisfy a closed set of lin¬ 
ear coupled equations as = C{{Om)) (where C de¬ 

notes a linear combination of On), then a two-time aver¬ 
age such as P{t)On{t + r) involving On will also evolve 
obeying the same set of equations, d{P(t)opt-\-T)) _ 

>C((P(t)Oni(t + r))). For example, the equation for the 
correlation function (cr+(t)cr“ (t-h r)) is obtained using 
the regression theorem in (24) as 

q 

+ {a+{t)a-{t + T)bg{t + T))P (29) 

Similar equations can be derived for other two-time 
operator averages to form a closed set, which can 
be numerically integrated starting from proper initial 
conditions given by quantities such as (cr+(ts)cr“(ts)), 
(cr+(p)cr“(p)6g(p))^, (cr+(p)a(p)), where p denotes 
the time when the system reaches steady state. The 
above steady-state expectation values can be evaluated 
using Bloch equations for these expectation values. The 
expectation values are also determined under the two- 
phonon correlation expansion approximation. As shown 
below (Sec. IIP.), these correlation functions will be used 
for calculating the WEA emission spectra and are accu¬ 
rate at low temperatures^. 
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Finally, it should be noted that a polaron shift Ap 
of the QD frequency (ujx) automatically appears in the 
spectrum calculated from the correlation expansion ap¬ 
proach and, for convenience, can again be absorbed into 
the polaron shifted QD frequency 


D. Emission spectrum 


1. Exact expression for the photon Green function spectrum 
for a generalized photon reservoir 


The incoherent emission spectrum of a QD coupled to 
a structured reservoir at a point detector at position rp 
is exactly given by— 

5'‘^(rD,w) = ((E+(rD,w))'''E+(rD,w)) , (30) 


where E'''(rD,w) = is the Laplace 

transform of the positive frequency component of the to¬ 
tal electric field operator E(rD,t), defined as 

E(rD,t) = E+(rD,t) + E“(rD,t) 

pOO 

/ (iw'[E+(rD, w',t) + E“(rD, 

Jo 

j j c^r'G(rD, r'; 

(31) 


/o 

+ H.C. 

In the frequency domain, 

E+ 


(rD,ct;) = i J duo' J dr'G{r-D^r']uj') 


X 


h 


Treo 




(32) 


Starting from the original Hamiltonian If in (JL) , the elec¬ 
tric field operator can be expressed in terms of the QD 
polarization using Laplace transform techniques, yield- 


4 Q 

mg— 


E(rD,w) = E°(rD,w) + G(rD,rd,w) • (33) 

where E^ denotes the free-held solution in the absence of 
a QD scatterer. The medium Green function G(rD, cj) 
includes all propagation effects^, including radiative 
losses due to structured reservoirs^; e.g., in the case of a 
planar geometry, this would also include out of plane ra¬ 
diative losses to photonic modes above the light line^. In 
the following treatment, we assume that the Green func¬ 
tion only accounts for radiative losses due to the struc¬ 
tured photonic reservoir. Any small additional out of 
plane radiative losses is accounted for by the phenomeno¬ 
logical Lindblad term 70 . Assuming the initial vacuum 
state of the photonic reservoir, the incoherent spectrum 
from Green function theory can be derived using f3Q) and 


(M), giving 


S‘^(rD,u]) 


G(rD,rd;uj) ■ d 

eo 


(( 7 +(w)cr (w)) 


= aprop(rD, rd; co) So(co), (34) 


where So = (a~^(uj)a~ (cu)) is the polarization spectrum 
and (aprop = 1 ^ ‘ C 1 (Rd, cj)!^ accounts for propaga¬ 

tion and hltering from the QD (at to the detector (at 
td). The polarization spectrum can be written as 


So(lo) = (a'^(oj)a (lo)) 

/ / POO \ t POO ^ 

^ \\J J dt2e^'^^‘^a~{t2) 

pOO pOO 

= dti dt2e^'^^^^~'^^^ {ti)a~ {t2)) . 

Jo Jo 


(35) 


Using a rotating frame at the exciton frequency (cj^) de¬ 
noting r = t 2 — ti and taking the limit of ti = t ^ oo, 
we derive the steady-state polarization spectrum. 


So{oj) = lim Re 

t^oo 


POO 

/ dr (cr+ {t + T)a- (t)} 

Jo 


(36) 


which can be calculated using the techniques developed 
in the previous sections. When calculated using the po- 
laronic approach (Secs. IIA. JIB.), a transformation is 
required to obtain So in the lab frame. If cTp denotes the 
QD lowering operator in the polaron frame, then a trans¬ 
formation of the correlation function ((cTp (t + r)crp (t))) 
from the polaron to the lab frame, 

(cr^(t -h r)crp (t)) ^ (cr^(t + r)B^{t -h r)B_{t)(T~ {t)) , 

(37) 


which produces a phase relaxation term that 

accounts for the phonon-induced pure dephasing of 
the QD polarization. The overall decay is thus 
clearly non-Markovian and this is a strong advan¬ 
tage of the polaronic approaches, i.e., even though 
they use a Born-Markov approximation for the equa¬ 
tions of motion, non-Markovian coupling to the phonon 
reservoir are captured through the polaron trans¬ 
form (indeed, the approach fully includes the IBM 
solution). The final polarization spectrum calcu¬ 
lated using the polaron approach is then So{oo) = 
limt^oo dr (cr^(t r)crp (t)) The 

expression for the emission spectrum (34) is exact and 
only limited by the approximations made in calculation 
of So from the above theories. Thus no phenomenologi¬ 
cal input-output formalism is required here as the Green 
function is already a solution for the scattering problem 
and the propagation of light in the medium. In order to 
distinguish S^ computed from the several different ap¬ 
proaches for obtaining the system dynamics, we identify 
the following spectra: and 
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when So is calculated using reservoir polaron ME, cQED 
polaron ME and the correlation expansion approach, re¬ 
spectively; the subscripts refer to a cavity or CROW 
medium as examples of the photonic bath, but note only 
the former approach can model arbitrary photon bath 
spectral functions (e.g., a non-Lorentzian lineshape). 

The propagator function (aprop(^^) filters the polariza¬ 
tion spectra So as the light propagated from the QD to 
the detector, and depends on the specific structure of 
the reservoir. Eor a single mode cavity with Lorentzian 
broadening, it is given by— 


'-^prop 


(w) 


fc(rD)®f*(rd)-d 

2 


eo 


— inuj 


(38) 


where fc(r) is the normalized transverse mode of the cav¬ 
ity. In the case of a high-Q cavity, and assuming k Uc, 
(38) reduces to a simple Lorentzian lineshape and the 
resulting Green function emission spectrum is given by 

'5Sv(rD,w) = ^ (cr+(w)cr“(w)) , 

[LU UJc) ~r [2J 

(39) 


where 


denotes the frequency- 


fc(rD)0f*(rd)-d ujI 
eo 2k 

independent prefactor. 

Eor a slow-light coupled-cavity waveguide (CROW)—, 
the frequency-dependent propagation factor is 


a; 


» = 


fcrow(l*D) C) f^j.Q^(rd) • d 


eo 


2 2 

Co’"' 


^/{uJ -UJu) 


(40) 


where uju,\ = <^u,i + ^eid cju,! marks the upper and 
lower mode-edge frequencies of the waveguide^!, de¬ 
notes damping and f^^Q^(r) is the normalized mode of the 
cavity containing the QD in the CROW—’—. The Green 
function emission spectrum from the waveguide is then 


can be approximately described by a single mode lower¬ 
ing operator a (strictly valid in high Q case) in the system 
Hamiltonian H (9). The damping of the cavity mode to 
the environment is described using the phenomenological 
decay tz. In this case, a calculation of the reservoir/cavity 
emitted spectrum can be obtained using 

= ^’(rd,rD)-x 
TT 

dT(at(t + T)a(t))e*(“--“’^ 

= F(rd,rD)-((a(w))ta(w)) , (42) 

TT 

where a{uj) is the Laplace transform of a(t), and 
^(i*d,i*D) is a frequency-independent geometrical factor 
to account for the propagation from the QD to the detec¬ 
tor position. The superscript CM denotes coupled mode 
formalism which treats the QD and cavity as coupled 
modes and the photonic reservoir in this case is assumed 
to be described a single cavity mode a. The structure 
(width) of the cavity determines its damping rate and 
any such reservoir effects arising due to finite cavity life¬ 
time is incorporated using phenomenological decay terms 
K for the cavity operator. An analytical expression for 
a{uj) can be derived using the Heisenberg equations for 
a. Starting from the coupled mode system Hamiltonian 
H (£), the Heisenberg equation for the cavity mode op¬ 
erator is 

da f ^ K\ _ ,, , 

— = - + ^ja-iga , (43) 

where the phenomenological cavity decay is included. 
Laplace transforming the a equation gives a(co’) = + 

/ where oP is the free field. Assuming an 

initial vacuum state for the cavity and substituting a{uj) 
in coupled mode spectra S^^{uj)^ we derive 

(44) 


lim Re 

t—yoo 


[L 


cG—res 
^crow 


(rD,w) = a'^J' 


VG - 


{a+{Lo)a (uj)). 


1 

COj*)(u) - UJu) 


2 

X 


(41) 


where = 


Lrow (rD)0C„^(rd)-d 

eo 


2 

. Eor a fixed detector 


location, the frequency-independent pre-factors 
containing the spatial mode functions do not influence 
the normalized emission spectrum. 


Coupled mode approaeh with eorrelation expansion and a 
simple eoupled eavity 

In the case where the structured reservoir consists of a 
Lorentzian cavity, the dynamics of the photon reservoir 


This functional form of the coupled-mode spectrum S'™ 
is similar to the Green function spectrum f39) . Hence 
for a fixed position of detector, the normalized spectrum 
derived from these two theories should match closely in 
the case of a Lorentzian cavity reservoir. This, however, 
happens only if all the approximations used in the cal¬ 
culations of (cr+(co’)cr“ (cj)) and (a^ (co’)a(co’)) are appro¬ 
priate and consistent. Eor example, the polaron cQED 
approach is derived using the Born-Markov approxima¬ 
tion with polaron transformed interaction terms. When 
applied beyond the range of the Born-Markov approxima¬ 
tion, the coupled-mode and Green function spectra may 
lead to different results (shown later). Hence, we will 
use this models to investigate the validity of the differ¬ 
ent theories in Secs. HIB and HI C. The two approaches 
where the coupled-mode spectrum can be computed are 
polaron cQED ME and correlation expansion approach. 
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To distinguish calculated using correlation expan¬ 

sion and polaron cQED ME, we label these as 
and respectively. We discuss the latter cou¬ 

pled mode spectra in details below. 


where Mi = M[ + iM{ and M 2 = MJ + iM^ (r and i su¬ 
perscript denotes real and imaginary parts, respectively); 
while the factors, Ni = and N 2 = 

7 cd 54 + + «A'), where 53 = 2 MJ - i{g' + 2M\) 

and Qi = 2M1 + i{g' - 2M{). 


3. Coupled mode approaeh with the eavity-QED polaron 
master equation 

The coupled mode spectrum (i.e., not using the exact 
Green function expression) in the polaron frame using 
the polaron cQED ME approach is given by 


^CM cQED^^^ ^ Firn)- ((ap(w))^ap(w)) = 
^ dr (al{t + r)ap {t )^ 


ii„, Re 

TT t^oc 


(45) 


where ap is the cavity lowering operator in the polaron 
frame. However the polaron transform does not influ¬ 
ence the cavity operator a = ap = e^ae“^, since P = 
cr+cr“ ^(^^“^< 7 )— • Thus the coupled mode spectrum 

remains unchanged when a transformation is made from 
the polaron to the lab frame (i.e.,((ap(cc;))’t'ap(a;)) = 
((a(a;))ta(a;))). 

A simple analytic expression for the coupled mode 
spectra can be derived in the WEA us¬ 

ing the simplified cQED polaron ME fl3b We start 
from Bloch equations fl9) in (a) and (cr“). The equa¬ 
tions for the two-time correlations {t)a~ {t -h r)) and 
(a’t'(t)a(t + r)) are the same as (19.) from the quantum 
regression theorem—. These equations are subsequently 
Laplace transformed— and solved to give the cavity- 
emitted spectrum 


oCM-cQED ^ ^ + 9c {a^O' )gg . . 

D{uj)C{uj)Pg,g, ’ ^ 


t - ret* 

where D{uj) = {uj — ^ ) -h , C{uj) = {uj — Acx' — 

+ peff 

^) -h and the steady-state expectation values 


(a’f'a)^^ and (aV are calculated using Bloch equa- 
tions for ( 0 * 0 )^^, (o'+cr-)^^, and {cr+a)^^. One 

obtains 


, . , P((A'2 - +2Re[ffi7V2]) 

CVss = -h-’ 

(47) 


t -\ _ (^i(^ + 70 ) - KN 2 ) + PN 2 


(aV ) = 

\ / ss 


(P + 7o)(A'2-|7,d|V(t^)2) ’ 


(48) 


where D = (A^^- |7cdf + )(Tg^(P+ 7 q) + )- 

2 Re[^i(Ai(P-^ 7 o)-/i:A 2 )], and TTotai = T^^ + T®^ is the 
total dephasing rate and A' = Acx' + A^^^ — A^^^ is 
the net QD-cavity detuning; also, gi = 2M[ — i{g' — 2M2) 


4- Cavity-emitted speetrum from a semielassieal linear 

suseeptibility theory 


In this section we discuss the cavity-emitted spectrum 
derived using the linear susceptibility theory—’—. The 
linear susceptibility approach is semielassieal and has 
been previously used for analyzing emission spectrum 
of a coupled QD-cavity system, specifically for investi¬ 
gating the role of phonons on the asymmetric vacuum 
Rabi-doublet^ and for connecting to experiments on off- 
resonant cavity feeding—’—; but its predictions and as¬ 
sumptions have never been compared against the quan¬ 
tum optical approaches starting with the full Hamilto¬ 
nian (Secs. HA. HB and H C) . In order to carry out such 
a comparison, here we briefly review the linear suscepti¬ 
bility approach. The theory assumes a time-dependent 
QD polarizability aspa,(t) = where 

T^c = 7 o + Id is the total dephasing rate of the bare 
QD and Px{0) is a constant prefactor. In order to incor¬ 
porate phonon effects to all orders, the time-dependent 
IBM phase (f is added to the Lorentzian decay of bare 
QD in the expression of Px{t). Next Px{t) is Eourier 
transformed to generate the frequency-dependent linear 
polarizability/susceptibility, x{^) = — 

ic^Ta,—ci;Eph(co’)), where Eph(co’) is the phonon self-energy. 
In the absence of a photonic reservoir, y describes the lin¬ 
ear response of a QD in presence of phonons to a weak 
pump. The imaginery part of y gives a measure of the 
linear absorption spectrum, which has the IBM spectral 
line-shape^’— (see Eig. 3., dashed line). Coupling to a 
structured photonic reservoir is included in this theory 
by assuming photon and phonon correlations to be inde¬ 
pendent. Thus the photonic reservoir only produces a SE 
enhancement (Purcell effect) and a photonic Lamb shift. 
Eor a single Lorentzian cavity (low or high Q), the QD 
susceptibility is 


x(w) = 


2lo' 


(K)^ - - iu)T^ - wl](w)) - 


-iujk) 


(49) 


Now assuming that the exciton is initially excited, the 
cavity-emitted spectrum is derived to be^’— 


;;^G—sus 


= KF(rd,rD)x 




(«)^ - - wS(w)) - CCE" 


iuju) 


( 50 ) 
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where F(rd,rD) denotes the frequency-independent ge¬ 
ometrical factor that depends on the QD and detector 
position—. The spectrum can be approximated as a pro¬ 
jection of QD susceptibility x(co’) (49.) by the Lorentzian 
cavity function (_38)22’— and is hence identified with the 
superscript G of the Green function projected spectrum. 
In the absence of phonon coupling (but accounting for the 
ZPL decay), the spectrum has the form of two coupled 
Lorentzians. With phonon coupling, the QD component 
becomes non-Lorentzian due to presence of a complex 
self-energy. 

Though applied in the case of a simple cavity here, the 
semiclassical linear susceptibility theory, in principle, can 
be generalized to arbitrary photonic reservoirs22. The 
theory also does not make any Markov approximations. 
However the self-energy S(c(;) derived using this approach 
depends on the QD dephasing rate T^,, which means that 
the model not self-consistent. The theory also cannot dis¬ 
tinguish between physically distinct pure dephasing and 
radiative damping processes, though (using a ME ap¬ 
proach) these are found to affect the linear spectra in an 
identical manner. The formalism also lacks a clear inco¬ 
herent pumping scheme and assumes an initially excited 
QD (i.e., (cr+cr“) (t = 0) = 0). This assumption however 
makes no difference in the calculated spectra as long as 
the WE A is valid (see Appendix A for details). Einally, 
the susceptibility theory cannot model any non-linear op¬ 
tical effects, which can be easily included in polaronic 
ME approaches (though the phonon-mediated scattering 
rates are no longer analytically solvable in general). 

At this point we have now introduced the various for¬ 
malisms that can be used for calculating the emission 
spectra of excited QDs coupled to structured photonic 
reservoirs in the presence of acoustic phonons. Before dis¬ 
cussing the results, we next introduce the specific phonon 
reservoir parameters which will be used in the subsequent 
calculations and we also discuss the computational pros 
and cons of the different numerical schemes. 


E. Phonon model and QD parameters 

So far in the theory we have not made any specific 
assumption about the structure of the phonon bath. 
Eor small epitaxial semiconductor QDs, the dominant 
phonon interaction arises due to deformation potential 
coupling to LA phonon modest. The expression for 

phonon coupling Xq for deformation potential coupling to 

^2 

LA phonons is Xq = 2 where pd is the 

material density, V is the material volume, Cg is the sound 
velocity, D = Dvai — Dcon is the difference in deforma¬ 
tion potential between the valence and conduction bands 
and cjp is the phonon cut-off frequency determined by the 
confinement length of the electron. The sum over the dis¬ 
crete phonon modes Xq can be converted into an integral 
using the density of states of phonons D(uj) = (- 2^)3 


as follows, 

poo 

/ dujD{uj)\{w) 

<1 -^0 

dwVVwJ Jpn{u}) (51) 

where the continuous phonon spectral function, Jpn(co’) = 

^2 

apUJ^e and phonon coupling o^p = 4 ^ 2 ^^!^^ — • The 

extra factor of ^/V appearing in f51i can be absorbed in 
the higher phonon correlations (e.g. {a~bq)), when this 
transformation is used in solving the correlation expan¬ 
sion equations (Sec. IIC) . Eor polaron and reservoir ME 
approaches, we use the continuum form of the phonon 
spectral function Jpn(<^) and with parameters for In As 
QDs, with ap = 0.06 ps^ and ujp = 1 meV, consistent 
with experiments^^. The precise value of these phonon 
parameters will not change any of the qualitative find¬ 
ings below, and they could be used, e.g., for fitting ex¬ 
periments on a particular QD. 


F. Numerical approaches 

In the above we have introduced three main theoret¬ 
ical approaches for analyzing the behavior of QDs cou¬ 
pled to structured photonic reservoirs in the presence of 
phonons, namely polaron MEs (cQED for a simple cav¬ 
ity and bath approaches), correlation expansion and the 
linear susceptibility technique. The linear susceptibil¬ 
ity technique is analytic and computationally the easi¬ 
est. Both the polaron MEs and the correlation expan¬ 
sion treatments involve simultaneous solutions for first- 
order coupled ordinary differential equations. The MEs 
solve coupled equations for the density matrix elements 
of the QD-photonic reservoir system and the correlation 
expansion solves coupled Bloch equations for QD-cavity- 
phonon correlations. The first set of equations (i.e., po¬ 
laronic ME) is much smaller, since the phonon compo¬ 
nent is eliminated using the Born-Markov approximation 
(Sec. JIB.). Eor the WEA applied here, the relevant ba¬ 
sis states for a QD-cavity system are |0I), |I0) and |00), 
where in \XC), the first (second) element denotes occu¬ 
pation of exciton (cavity). Eor this purpose, we use the 
Quantum Optics (QO) Toolbox— to numerically solve 
the ME, which includes the states |02), |20), |II), for the 
truncation of the ME^. The total density matrix is then 
expressed in a basis of 6 states. The two-time correla¬ 
tion function O^t + r)0{t) for any operator O can be 
obtained by solving coupled first-order differential equa¬ 
tions for 16 density matrix elements. The 16 elements 
consist of 6 population elements (diagonal) and 10 po¬ 
larization elements (first off-diagonal)—. The forward 
matrix propagating the system in time is thus 16 x 16 
in size. When the polaron reservoir treatment (Sec. II A) 
is used, after elimination of the photonic reservoir, the 
number of basis states are those of the QD (|I) and |0)) 
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and the forward matrix is just 4 x 4 in size. In compari¬ 
son, the non-Markovian correlation expansion treatment 
scales as N‘^ in number of correlation elements and N‘^ 
X in forward matrix size, where N is the total num¬ 
ber of phonon modes used in discretizing the relevant 
regions of the phonon bands (±5 meV for the current 
QDs, Sec. II Eb Finer discretization reduces the numer¬ 
ical noise and for N = 400 (typical numbers), the for¬ 
ward matrix size is 160000 x 160000. Thus on a regular 
work station, the correlation expansion calculations can 
take a few days to get the cavity-emitted spectrum com¬ 
pared to polaron ME calculations which take only a few 
seconds and are this much easier to work with from a 
numerical perspective. Moreover, the ME approach can 
easily include other pump field and nonlinear processes, 
including multiphoton effects. The correlation expansion 
is however useful, if one works at low temperature, to 
justify and test the limits of the polaron ME approaches. 


III. RESULTS 

A. Emitted spectrum from an uncoupled QD 

We begin this results section with the spectrum of a 
simple QD exciton coupled to phonons, which is decou¬ 
pled from any structured photonic reservoirs. We ex¬ 
pect the spectrum to represent the IBM lineshape due to 
the presence of phonons with a broadened ZPL, which 
is a well known result and a useful starting point with 
which to understand the results when a photon bath is 
included. In this limit, and at low temperatures, all the 
models produce qualitatively identical spectra so we only 
use the polaron ME model here to analyze the results in 
this section. The QD is assumed to be weakly excited 
by an incoherent pump (P) and has a background de¬ 
cay rate 70 , and pure dephasing rate, jd- In the polaron 
frame (Sec. IIBb the QD emission spectrum (36) is given 
by 

pOO 

Sq{uj) = lim Re[ / dr (dp (t + r)crb (Q) 

t^OO Jq 

(52) 

which is a simple Lorentzian with a total ZPL linewidth 
7o + 7^5 as shown in Fig. 3_ (dash-dotted line, without 
phonon coupling). The transformation from the polaron 
frame to the lab frame produces a phase term (37) 
and the final QD emission spectrum (i.e., in the lab 
frame) is 

pOO 

So{uj)(x lim Re[ / dr (^ap {tr)ap {t)) 

t^oo Jq 

(53) 

The resulting emission spectrum (thick orange solid line) 
is plotted in Fig. 3_ and shows the appearance of the 
phonon sidebands arising due to the IBM phase func¬ 
tion (j). The asymmetry of the sidebands arise due to 



FIG. 3: (Color online). Normalized QD emission (polarization) 
spectrum (thick orange solid) and linear absorption spectrum (thick 
black dashed) at T = 4 K for a single uncoupled QD. The red dash- 
dotted line denotes the polarization spectrum Sq in the polaron 
frame. The ZPL parameters are 70 = 5 /reV, rd = ^ fieV. 


the fact that phonon emission is more probable than ab¬ 
sorption at low temperatures, so the phonon emission is 
more probable on the lower energy side of the ZPL. The 
Lorentzian ZPL corresponds to the polarization spectrum 
Sq in the polaron frame. 

Similar to Sec. IIP4. the QD susceptibility function- 
can be defined as 


pOO 

x{od) oc i lim / dr (cTp {t + r)ap (t )) 

t^OO Jq 




(54) 


and the linear absorption spectrum (Im(x)) is plotted 
in Fig. 3_ (thick black dashed line). The linear absorp¬ 
tion spectrum is simply a reflection of the emitted spec¬ 
trum about the QD ZPL^. The phonon sidebands are 
now more enhanced to the right, since phonon abso- 
prtion is more probably on the higher energy side of 
the ZPL. To explain this feature, consider a weak co¬ 
herent laser that excited the QD, i.e., through i^pump = 
rjx{o'~ where ujl and r]x are the drive 
strength and the Rabi frequency of the laser, respectively. 
A photon at higher frequency {ujl > can excite the 
QD more easily by phonon emission process at low tem¬ 
perature. Since the reverse process {ujl < ojx) requires 
phonon absorption, which is less probable at low tem¬ 
peratures, then absorption is stronger to the right (high 
energy side) of the ZPL. The linear absorption spec¬ 
tra will be used later in an attempt to understand off- 
resonant cavity feeding using linear susceptibility theory 
(See, me 2 b 

We next consider coupling to structured reservoirs and 
begin with the case of a resonant high-Q cavity, i.e., in 
the strong coupling regime. 
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FIG. 4: (Color online). On-resonance strong coupling 
regime. Cavity emission spectra at T = 4 K for a strongly cou¬ 
pled QD-cavity system at resonance = ujc)- The plots compare 
coupled mode spectra calculated using different approaches. The 
thick orange (light solid) line in (a), (b) and (c) denotes the coupled 
mode spectrum, calculated using correlation expansion ap¬ 

proach. The magenta (thin solid) and the black (dark) dashed line 
plots the coupled mode spectrum, cQED using po- 

laron cQED approach, in the absence and presence of phonon cou¬ 
pling in (a) and (b), respectively. Panel (c) compares the results 
with linear susceptibility theory and the black (dark) dashed line 
denotes the semiclassical Green function spectrum calcu¬ 

lated using a linear susceptibility approach. The main parameters 
are g = 100 /reV, = 65 p,eV, 70 = 5 p-eV, and 7 ^ = 55 p,eV. 


B. On-resonance strong coupling regime in a 
high-Q cavity 

As an example of strong coupling between a QD and 
a photonic reservoir, we consider the case a QD strongly 
coupled to a high-Q cavity, where the cavity is described 
by the cavity mode operator a. This is studied in Fig. T 
As is well known, a strongly coupled QD-cavity sys¬ 
tem undergoes vacuum Rabi oscillations, when a sin¬ 
gle quanta of energy is coherently exchanged between 
the QD and the cavity-’-. For a simple two-level atom, 
coupled to a symmetric cavity without phonon effects, 
the emission spectrum shows the two hybridized polari- 
ton states of equal intensity (magenta, thin solid) line. 
Fig. 4.(a), separated in frequency by twice the QD-cavity 


coupling constant g. Such a coherent transfer of en¬ 
ergy between light and matter is important for building 
quantum light-matter interfaces, which could be used, 
e.g., for long distance quantum communication with pho- 
tons2. Phonon interactions however affects the coherence 
of the system. In the presence of phonon coupling, the 
polaritons appear with different intensities and reduced 
vacuum Rabi-splitting (^ 2 {B) orange (thick solid) 
line. Fig. ^ (a))—’—. Such phonon-dressed strong cou¬ 
pling can be well explained using the polaron cQED the¬ 
ory (Sec. JIB.), the correlation expansion (Sec. IIC) or 
the linear susceptibility (Sec. IIP 1) approaches, and we 
compare these different spectra in Fig. T In Fig. 4.(b), we 
compare the coupled mode spectra derived using the cor¬ 
relation expansion solid line) and cQED polaron 

ME technique dark dashed line) and the re¬ 

sultant spectra show very good agreement. The QD (cj^) 
and cavity (cjc) are assumed to be at a detuning of 0 
meV (ZPL) and the frequency of the QD, includes 
the polaron shift. The parameters used for the simu¬ 
lations are close to typical experimental numbers, which 
show vacuum Rabi splitting of 2g = 200 //eV— {2g' = 183 
/ieV) in Eig. T As discussed earlier, the polaron cQED 
approach is derived using a Born-Markov approximation 
for the polaron-transformed interaction terms and is valid 
as long as the system dynamics (i.e., the vacuum Rabi 
period tt/^ ~ 20 ps (Eig. 4.)) is longer than phonon re¬ 
laxation time (~ 3 — 5 ps. Sec. HE) . The correlation 
expansion approach however uses no such approxima¬ 
tions and is expected to work as long as the temperature 
is sufficiently low (e.g., less than 20 K or so), though 
it is numerically more cumbersome as discussed earlier 
(Sec. HE) . The close agreement (Eig. 4.(b)) of the polaron 
cQED results (dark dashed line) with the more rigorous 
non-Markovian correlation expansion results (solid line) 
demonstrates the validity of polaron cQED ME technique 
which is also significantly less numerically demanding. 

In Eig. 4.(c), we also compare the coupled mode 
spectrum obtained from correlation expansion approach 
(solid line) and the semiclassical Green function 
spectrum obtained using linear susceptibility approach 
(dashed line) and they show very good agreement. 
As described in Sec. IIP4. the semiclassical linear sus¬ 
ceptibility approach does not make any Markov approx¬ 
imation and evidently works well when the on-resonance 
condition {/S^cx' = 0)—’— is satisfied. As we will show 
in the next section, however, this semiclassical approach 
fails to get the correct spectra in the off-resonant cavity 
feeding regime {/S^cx' ^ ^)- 

As seen above, one of the most striking effects from 
phonon coupling at the on-resonance {/S^cx' = 0) condi¬ 
tion is the asymmetric vacuum Rabi doublet22’AI. The 
exact form of the ME (12) however does not bring into 
light the physical processes responsible for this phonon 
induced assymetry. The same problem lies with the cor¬ 
relation expansion approach which mathematically re¬ 
duces to a set of coupled first-order differential equations. 
The full polaron cQED ME (12) can however be sim- 
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FIG. 5: (Color online). Normalized emission spectra at T = 4 
K for a strongly coupled QD-cavity system at resonance. The plots 
compare the Green function spectra (solid line) and the cou¬ 
pled mode spectra (dashed line) calculated using correlation 

expansion (top panels (a, c)) and polaron cQED ME approaches 
(bottom panels (b, d)). The QD-cavity coupling g is 50 fieV for (a, 
b) and 100 p,eV for (c, d). The other parameters are k, = 65 
70 = 5 fieV, = 55 fieY. 


plified to an analytical form (13) . as long as the WE A 
is valid. This is true for the current situation as we 
are only dealing with the linear spectrum (as shown in 
Fig. (4.)). The analytical cQED ME form (13) consists 
of incoherent feeding terms (T^^^ /cr+a^^ cross-dephasing 
term 7cd, and additional terms denoted by Mi, M 2 . For a 
weakly excited system at resonance, the incoherent feed¬ 
ing terms and the cross-dephasing term 7cd have negligi¬ 
ble effects on the spectrum. Thus the term solely respon¬ 
sible for the asymmetric doublet is then M 2 (as Mi = 0 
when Acx' = 0). Such a term also causes an asymmetry 
in the ME approaches that describing high-held Mollow 
triplets^ when they sample an asymmetric bath. This 
term then acts as an effective complex coupling ig' — M 2 
between the resonant QD and cavity (19) . 

The polaron cQED ME approach is strictly valid when 
the Markov approximation is satished and, as was shown 
above, the coupled mode spectra calculated using the 
polaron cQED approach matches the more rigorous cor¬ 
relation expansion results closely. It was also mentioned 
earlier in Sec IIP2. that the Green function spectrum 
*5^^ and the coupled mode spectrum S'™ derived us¬ 
ing the same theory should match, as long all underlying 
approximations are valid. In order to test this, we ap¬ 
ply non-Markovian correlation expansion and cQED ME 
techniques to calculate low temperature (T = 4 K) vac¬ 
uum Rabi spectrum for parameters which do (Fig. A (a, 
b))- and do not satisfy (Fig. A (c, d))— the polaronic 
Born-Markov approximation. The parameters used for 
both simulations are close to typical experimental num¬ 
bers, which show vacuum Rabi splitting of 2g = 100 
/reV- {2g' = 91.5 geV) for Fig. A (a, b) and 2g = 200 


geV— {2g' = 183 /reV) for Fig. A (c, d). As expected, 
the Green function spectrum (thick solid line) and 

the coupled mode spectrum (dark dashed line), 

calculated using correlation expansion match exactly in 
both cases (Fig. A (a, c)). For Fig. A(b), the vacuum 
Rabi period ir/g ~ 40 ps is much slower than phonon 
damping (4 ps), which makes the Markov approxima¬ 
tion valid. Thus the Green function spectrum 
(solid line) and the coupled mode spectrum 
(dark-dashed) from the polaron cQED approach match 
closely. However for Fig. A(d), the vacuum Rabi pe¬ 
riod ir/g ~ 20 ps, and is now closer to the phonon 
relaxation time, the Born-Markov approximation may 
fail. Hence the coupled-mode spectrum (dark 

dashed line) now differs strongly from the Green func¬ 
tion spectrum (thick solid line), which shows 

an appearance of an unphysical third peak between the 
two polariton peaks (Fig. A(d)). This happens due to 
the break-down of the Born-Markov approximation, and 
the polaron cQED approach over-estimates the phonon 
dressed emission from the QD, at frequencies between the 
polariton peaks. Hence, when using the polaron cQED 
theory, the physically correct spectra is obtained using 
the coupled mode approach. As shown above, it is bet¬ 
ter to use the coupled mode approach of the polaron 
cQED ME since the cavity emitted spectrum is consistent 
with the input-output theory, while the Green function 
approach uses the exact input-output relations and can 
bring back extra features that are not contained in the 
Born-Markov approximation. 

It should be noted that a rather large value of pure 
dephasing (= 55 geY) is used for these spectra cal¬ 
culations. Such a large value is chosen to achieve faster 
numerical convergence for the correlation expansion cal¬ 
culations, which are numerically intensive. This value of 
Yd is also used for the off-resonant cavity feeding calcu¬ 
lation using correlation expansion, in the next section. 
The qualitative effects for different values of jd remain 
the same. 


C. OfF-resonant cavity feeding 

1. Comparison of three theoretieal approaehes 

In this section we consider the case of phonon-mediated 
cavity feeding from an off-resonant QD22. Because of the 
off-resonant condition {Acx' > 1 nieV), an approximation 
of weak-coupling between the QD and cavity, is not very 
restrictive. This allows us to use the polaron reservoir 
theory (Sec. HA) to calculate the spectrum along with 
the polaron cQED ME approach, correlation expansion 
and semiclassical linear susceptibility theories. As shown 
earlier, the reservoir theory accounts for dynamical inter¬ 
play between photon and phonon baths, which is absent 
(or restricted to certain regimes) in the other theories^. 

Off-resonant feeding is investigated in the case of a 
high-Q (Fig. 6.) and low-Q cavity (Fig. 7), located ap- 
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proximately - 2 meV away from the QD. In the pres¬ 
ence of phonon coupling, the intensity at the off-resonant 
cavity is strongly enhanced (Fig. 6 _(a), solid line). Such 
strong feeding is not expected without phonon coupling, 
when only pure dephasing of the QD is accounted for 
(Fig. 6 .(a), black dashed line). The feeding is much 
smaller in the low-Q case (Fig. i.(a), solid line)) and we 
magnify the region near the origin (Fig. 2.(a), dashed line) 
to better highlight the cavity feeding effect. The sub¬ 
sequent graphs in the low-Q case (Figs. i.(b, c)) focus 
near the origin to investigate off-resonant feeding. We 
compare spectra calculated using the different theoreti¬ 
cal approaches in Fig. 6 . and Fig. 7., where Note the Green 
function spectrum, and the coupled mode spectra 

match exactly when calculated using the non- 
Markovian correlation expansion approach. Hence we 
use the coupled mode spectrum, (solid orange 

(thick light) line. Fig. 6 . and Fig. i), to compare against 
the other approaches. Figures 6 . and 7_ (b) compares 
the correlation expansion calculations with the polaron 
cQED ME approach (dashed line). Eollowing the discus¬ 
sion in Sec. IllB. the coupled mode spectrum 
from the polaron cQED approach is chosen for compar¬ 
ison. The spectra compare closely in the high Q cavity 
(Eig. 6 . (b)) and show large disagreement in the low Q 
cavity (Eig. 7. (b)). Eigures 6 . and 7. (c) compares the 
correlation expansion spectrum with the Green function 
spectrum derived using the polaron reservoir approach 
black (dark) dashed line) and the results show 
excellent agreement in both high and low Q cases. Thus 
the coupled mode spectra of the polaron cQED approach 
is not well suited for low Q cavities. 

The emission spectrum (solid orange line) resembles 
a sum of two coupled Lorentzians in the high Q cavity 
(Fig. £.(c)). In the low Q cavity however, the emission 
spectra has a non-Lorentzian shape (Fig. Z(c), solid or¬ 
ange line) which resembles the IBM spectral lineshape. 
Physically, for low Q cavities the emission spectrum 
should resemble the non-Lorentzian IBM spectrum as 
demonstrated experimentally by Weiler et alJ^, where 
low Q cavities are used to measure vertical emission from 
QDs. The Green function spectrum from the reservoir 
theory correctly reproduces the spectrum (black 

(dark) dashed line Fig. Z(c)) in the high and low Q 
case and the reason for this can be explained by under¬ 
standing the way in which it is calculated. The calcu¬ 
lation of the Green function spectrum using the reser¬ 
voir approach uses the polarization spectra So (53) fol¬ 
lowed by subsequent projection using the cavity prop¬ 
agator aprop (-SS.), which is an exact input-output for¬ 
malism. The SE rate of the QD is given by 7 — us¬ 
ing the photon reservoir approach (^). The exciton po¬ 
larization spectra in the polaron frame f52) . Sq (uu) = 
limt^oo (<^p (^ + (^))] thus a Lorentzian 

(dark dash-dotted line, Eig. ^) with spectral width de¬ 
termined by total dephasing 7 + 70 + 7 d + where the 
latter is negligible here. This Lorentzian (dash-dotted 
line, Eig. ^) is the ZPL of the QD. The non-Lorentzian 
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FIG. 6: (Color online). Off-resonance cavity feeding 
regime. Normalized emission spectra at T = 4K for a QD weakly 
coupled to a high-Q cavity (Q = 8000 at =1440 meV), un¬ 

der off-resonance condition. The thick solid orange (light) line in 
(a), (b) and (c) denotes the coupled mode spectrum, cal¬ 

culated using correlation expansion approach. The black (dark) 
dashed line plots the coupled mode spectra, in the ab¬ 

sence (presence) of phonon coupling in (a) ((b)) using the polaron 
cQED ME approach. Panel (c) compares the correlation expansion 
calculations (thick solid line) against the Green’s function spectrum 
calculated using polaron reservoir approach black (dark) 

dashed) and linear susceptibility (S'^ 7^^^5 magenta (light) dashed) 
approaches, respectively. The parameters are g = 100 /reV, k, = 
180 /reV, 70 = 5 /i-eV, = 55 yaeV, A^a,/ = 2 meV. 


sidebands (thick light solid line. Fig. ^) of the QD emis¬ 
sion spectrum in the lab frame So (53) due to acous¬ 
tic phonon interaction— arise as the phonon correlation 
function is analytically included in the QD corre¬ 

lation (37) when a transformation is made from polaron 
{So) to lab frame {So) (see discussion in Sec. Ill A) . This 
gives So the characteristic non-Lorentzian lineshape of 
the IBM (Fig. Z, light solid line) and is plotted for the 
low Q case in Fig. Z{c) (blue dashed-dotted line)—. The 
final Green function spectrum (Fig. i.(c), black 

(dark) dashed line) is obtained by multiplying the polar¬ 
ization spectrum So (Fig. i.(c), blue dashed-dotted line), 
with the cavity propagator <aprop (38), which amplifies 
the phonon sidebands. In the high Q case, the cavity 
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FIG. 7: (Color online). Normalized emission spectra at T = 
4 K for a QD weakly coupled to a low-Q off-resonant cavity (Q 
= 600 at a;c/27r =1440 meV). The thick solid orange (light) line 
in (a), (b) and (c) denotes the coupled mode spectrum, 
calculated using correlation expansion approach. The black (dark) 
dashed line in (a) shows phonon induced feeding and is obtained 
by magnifying (x5) the spectrum in this region. Figures (b, c) 
zoom in on the region near origin, to magnify the phonon effects. 
Plot (b) compares correlation expansion calculation (solid orange 
line) with coupled-mode spectra cQED (dark) dashed) 

calculated using polaron cQED ME approach. Plot (c) compares 
the correlation expansion calculations (solid orange line) against 
the Green function spectrum calculated using polaron reservoir ap¬ 
proach black dashed) and linear susceptibility 

magenta (light) dashed) approaches respectively. The blue dash- 
dotted line shows the polarization spectra Sq calculated using the 
polaron reservoir approach which resembles the IBM spectra. The 
parameters are g = 100 /ieV, k = 2.4 meV, 70 = 5 /ieV, jd = 55 
/reV, = 2 meV. 


propagator amplifies a narrow region of the phonon side¬ 
band around oOc and the Green function spectrum from 
the reservoir approach resembles a sum of two 

Lorentzians at QD and cavity. This is however not the 
situation in the low Q cavity (Fig. 2.(c)). For a low- 
Q cavity, the position of the second peak due to cavity 
amplification strongly depends on the phonon sideband 
spectrum and may appear at a position different from 
its original location at cjc ~ meV— . The process 
also makes the resultant spectra (black (dark) 

dashed line) in Fig. 2.(c) resemble the IBM spectrum {Sq^ 


blue dash-dotted, Fig. 7.(c)) more closely. Thus the pro¬ 
jected spectra from the photon reservoir theory matches 
the physically correct results closely and provides a huge 
computational advantage over the more numerically de¬ 
manding correlation expansion technique. 

The coupled mode spectrum from the po¬ 

laron cQED approach resembles a sum of two Lorentzian 
modes with the QD and the cavity peaks appearing at 
their original positions = 0 and ujc ^ 2 meV, respec¬ 
tively (Fig. 6 .(b) and Fig. 2.(b), black dashed line). This 
is close to the correct spectral shape in the high Q case 
(solid orange line. Fig. £.(b)) but deviates strongly from it 
in the low Q case. In the low Q case, the non-Markovian 
dynamics is not correctly accounted in polaron cQED 
coupled-mode spectra oc (a’*'(cc;)a(cc;)), which 

always resembles a sum of two Lorentzians. This non- 
Markovian dynamical effect is captured by the correla¬ 
tion expansion calculation as well as the Green function 
spectrum from the reservoir approach. In the reservoir 
approach, this dynamics is accounted in the QD polar¬ 
ization spectra Sq, which is subsequently multiplied by 
the Lorentzian cavity projector to calculate 

A slight disagreement between the projected reservoir 
spectrum (black dashed) and coupled mode cor¬ 
relation expansion spectrum (orange solid) can 

be observed in the low Q case (Eig. 7.(c)) as compared 
to the high Q case (Eig. £.(c)). This difference arises, 
because the reservoir theory accounts for the dynamical 
interplay between the photon and phonon bath in de¬ 
termining the SE rate 7 (^); 7 in turn determines the 
spectral-width of the ZPL in the QD polarization spec¬ 
tra So (53) . The structure of the photonic reservoir is 
more correctly incorporated in the reservoir approach as 
compared to correlation expansion technique, which in¬ 
corporates the photon bath at the level of a system op¬ 
erator a and accounts for its structure / damping us¬ 
ing phenomenological decay term tz. In the latter case 
7 = 7 = (21). Such an ap- 

proximation is strictly valid for a high Q cavity, where 
the spectral width of the photon bath is negligible com¬ 
pared the phonon reservoir band-width ( 8-10 meV) and 
breaks down (7 7 ^ 7 p) when both reservoirs have com¬ 
parable spectral widths. This is the case for the low Q 
cavity. The difference is however very small here because 
the net dephasing rate 7 + 70 + 7 d + determining the 
spectral width of the Lorentzian ZPL of So is dominated 
by the pure dephasing ( 7 ^ = 55 fieV) in the current case, 
which is much larger than 7 1-4 jgeV. 

A comparison between the spectra from the correlation 
expansion approach (light solid), and linear sus¬ 
ceptibility theory (magenta (thin light) dashed) is 

also shown in Eig. £.(c) and 7.(c), where the linear suscep¬ 
tibility theory shows very little feeding to the off-resonant 
cavity. This discrepancy with linear susceptibility theory 
in the context of off-resonant feeding was mentioned ear¬ 
lier in Sec. IIIB and is explained in detail in Sec. Ill C 2 
below. 
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FIG. 8 : (Color online). Off-resonance cavity feeding regime: polaron ME versus linear susceptibility model. Normalized 
Green function spectra for a QD-cavity system without phonons (black dashed) and with phonons (T = 4 K) (solid orange) calculated from 
linear susceptibility method (fop panels) and the polaron reservoir method, (bottom panels). The QD-cavity detunings 

are 1, 2, 3, 4 meV and -1, -2, -3, -4 meV (top to bottom) in left and right panels respectively. Individual spectra are normalized by peak 
ZPL intensity and data in bottom-left panel is multiplied by 1.5 for better visibility. The main parameters are g = 100 /reV, k = 65 /reV, 
70 = 5 /reV, 7 ^ = 55 /reV. 


2. Explanation of the failure of the susceptibility theory for 
explaining the off-resonant cavity feeding 


the linear suspecpibility approach to obtaining the spec- 


In this section we compare the Green function spec¬ 
trum calculated using the polaron reservoir ME (^) and 
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trum (50) . for the case of off-resonant cavity feeding 
and explain the large differences between the two re¬ 
sults shown previously in Fig. 6. and 7_{c). As shown in 
Ref. 33^ the cavity-emitted spectrum (i-^- -SO.) is 

obtained by multiplying the linear susceptibility y with 
a Lorentzian cavity projector (38) . The Im(x) (linear 
absorption spectra) plotted in Fig.^ (dark dashed line) 
shows that the phonon induced enhancements (asymme¬ 
tries) are stronger to the right, than to the left of the 
ZPL. The same holds true for the Re(x) (not shown). 
Thus when multiplied by a Lorentzian projector for the 
spectrum f5Q) , cavity feeding is stronger to the right than 
to the left^’— . This is seen from the graphs in the top 
panel of Fig. where normalized spectra (light solid 
orange line) are plotted for different detunings in the 
presence of phonons. The dark dashed lines represent 
the corresponding normalized spectra in the absence of 
phonons. Individual spectra are normalized with respect 
to the peak ZPL intensity. Left (right) panels plots spec¬ 
tra for positive (negative) QD-cavity detunings (Ax'c)- 
Cavity feeding can be estimated from the ratio of peak 
heights as both Lorentzians have comparable linewidths 
(Tx = 60 /reV and k = 65 /reV) and it increases in pres¬ 
ence of phonons. The feeding is asymmetric and stronger 
when the cavity is to the right. These results are how¬ 
ever in complete contrast with the polaron reservoir ME 
calculations (bottom panels, solid orange lines), which 
shows a stronger cavity enhancement to the left of the 
QD ZPL. The dark dashed lines once again plot the nor¬ 
malized spectra without phonons. As explained earlier 
(Sec. IIP lb the Green function spectra, (34j is 

obtained by multiplying the polarization spectrum Sq of 
the QD (Fig. solid line), with the Lorentzian cavity 
projector (38) . The QD polarization spectrum (Fig. 
solid line) shows strong phonon bands to the left of the 
QD ZPL and hence the cavity is fed more strongly to the 
left at low temperatures. 

We have so far provided a mathematical explanation 
for this difference in cavity feeding. The results from 
the polaron reservoir ME can however be justified to be 
physically correct. At low temperatures phonon emis¬ 
sion is more probable than phonon absorption. Thus 
if a cavity has lower energy than an excited QD, the 
cavity can be excited more easily in a two-step quan¬ 
tum process where a cavity photon is created along with 
a phonon emission. When Uc > exciting the cav¬ 
ity will require phonon absorption which is less feasible 
at low temperatures. Thus the polaron reservoir ME 
predicts the physically correct solution. The linear sus¬ 
ceptibility theory misses this two step quantum process 
which is correctly incorporated in the SE rate f21) cal¬ 
culated using polaron reservoir ME approaches, through 
the Lindblad QD/cavity scattering rates . The 

linear susceptibility technique is however very successful 
in explaining the on-resonance {Acx' = 0) phonon in¬ 
duced assymetric vaccum Rabi doublet . This is because 
at resonance, these two-step quantum feeding processes 
barely contribute to the emission spectra (see Sec. IIIBb 


The technique however fails for far detuned condition 
{\Acx'\ ^ 0), where these two-step quantum processes 
become important. 


D. Emission spectra from a slow-light 
coupled-cavity waveguide 

The reservoir theory can be applied to general struc¬ 
tured reservoirs, beyond simple Lorentzian cavities. In 
photonic structures with a non-Lorentzian LDOS pro¬ 
file, the dynamical interplay between photon and phonon 
reservoir can manifest itself more strongly, as has been 
demonstrated by photoluminescence intensity studies in 
coupled-cavity waveguides^. In this section we calcu- 
ate the emission spectra of a QD coupled to a 

photonic crystal waveguide, using the reservoir theory. 
Photonic crystal waveguides are important for slow light 
applications^”— and on chip single photon emission—”—. 
The specific case of a photonic crystal coupled-cavity 
waveguide (CROW)^ is considered here for which a sim¬ 
ple expression of the medium electric field Green func¬ 
tion, G can be derived analytically, using a tight-binding 
approach^”—. The Green function is in turn used to 
determine the LDOS and the projector ap (40) of the 
wave-guide. The Purcell factor (PE) and the projector is 
plotted in Eig. £(a) and (b), respectively. The PE (which 
depends on the projected LDOS) is given by 7 / 7 ^,, where 
7 = 2 Re[Jph(T)]dr— is the frequency-dependent SE 
rate of the QD in the waveguide without phonon cou¬ 
pling, and Jph(^) is the waveguide bath relaxation func¬ 
tion (Sec. IIAj and 7 ^ = d‘^^/euJ^/{dirheoc^) is the de¬ 
cay rate in the background slab material. As shown in 
Eig. £(b), the projector can be approximated as a sum 
of two sharp Lorentzians located at the upper (uju) and 
lower {uji) mode-edges of the waveguide. The parameters 
used correspond to a GROW structure built by coupling 
individual cavities formed by local width modulation of a 
line defect waveguide^, on a photonic crystal slab, which 
has a period of d = 420 nm — . The individual cavities 
have a mode volume of Wff = 0.175 x 10“^^ m^ and the 
GROW period is 5d^’^. 

The polarization spectrum from a bare QD resembles 
the IBM spectruml2’— (Eig. solid line). In the weak- 
coupling limit, the spectrum of a QD coupled to waveg¬ 
uide (Eig. 10 ) is obtained by multiplying the IBM spec¬ 
trum with the projector (41) . Eignres 10 (a) and (b) plot 
the waveguide spectrum (normalized units) at T 

= 4 and 40 K respectively, for a QD located at the band- 
center, ujQ {= As shown in Eig. 10 (b). the projector 
amplifies the phonon side-bands and leads to the appear¬ 
ance of three distinct emission peaks, which corresponds 
to the two waveguide mode-edges and the QD ZPL. The 
peak intensity of the three peaks can be denoted as Iq = 
57---K) (QD-ZPL peak), Uq = (upper 

waveguide band (UB) at cj = uju) and Lq = 

(lower waveguide band LB at uj = uJi). The emission 
at the lower mode-edge (Lq) is stronger than the upper 
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FIG. 9: (Color online). Slow-light CROW Purcell fac¬ 
tor and propagator. (a) Purcell factor (PF) and (b) pro¬ 
jector function ap for a coupled-cavity photonic crystal waveg¬ 
uide (CROW). The parameters used for waveguide calculations 
are taken from Ref. ^ (see text) and QD-waveguide coupling 

g = (^ 2heQ^V ^ = 85 geY (with = 50 Debye), where ouq marks 
the waveguide band-center. 
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mode-edge (t/o), primarily because of the asymmetry of 
the phonon sidebands due to a stronger phonon emission 
probability at low temperaturesl^. The phonon-assisted 
“feeding” of the waveguide mode-edges can be varied by 
tuning the QD in frequency as shown by Fig. lQ(c). which 
plots (at T = 40 K) for a QD located 1 meV to 

the right of the band-center (ujq). As seen from the nor¬ 
malized spectra of Fig. 10_ (a) and (b), mode-edge feeding 
due to phonons become stronger at higher temperatures. 
Hence we use a higher temperature (T = 40 K) for further 
investigating phonon effects below. For the moment, any 
temperature dependence of pure dephasing is neglected 
(Fig. 10), but we will include this effect below (Fig 13(b)). 

As also observed from Fig. 10., the spectral width of the 
waveguide (8 meV) is comparable to the phonon reser¬ 
voir. The dynamical interplay between the two reser¬ 
voirs is expected to influence the emission spectra in this 
case. As shown elsewhere^’—, this interaction leads to 
a strong enhancement of SE rates outside the waveguide 
band (also Fig. 12.(a)), where SE otherwise is strongly 
suppressed (Eig. £(a)). Eor the case of the waveguide 
emission spectrum, the effect of this interplay can be 
observed by tuning the QD close to the upper mode- 
edge (Eig. 11 ) . Since the lower mode-edge is located at 
a frequency larger than the spectral width of the phonon 
side-bands, its emission intensity is negligible and hence 


FIG. 10: (Golor online). Normalized emission spectra 
for a QD weakly coupled to a waveguide, (a) and (b) plot the 
spectrum at T = 4 and 40 K, respectively, when = cjq, where 
(jjQ marks the waveguide band-center and (c) is the spectrum at T 
= 40 K and uo'^ — ujq = 1 meV. The UB and LB mark the upper 
and lower waveguide bands. The parameters are = 85 /i-eV, 70 = 
7d = 1 geV. 


not included in the current figure. In Fig. JT (a), the 
QD (at 0 meV) is tuned 1 meV outside the upper mode- 
edge. The upper mode-edge appears clearly on the spec¬ 
trum due to the presence of the phonon sidebands. The 
solid line plots the emission spectrum, when the SE rate 
(7) of the QD is estimated by taking the effects phonon 
bath into account (A). The dashed line plots the spec¬ 
tra, for the case when phonons do not influence the SE 
rate (7 = 2 Re[Jph(r)](ir)H. Due to the phonon- 
mediated SE enhancement outside the wave-guide mode 
edge, the solid line shows a stronger feeding of the waveg¬ 
uide compared to dashed line. The reverse effect happens 
in Eig. JT (b) when the QD is tuned inside the waveguide 
band by about 0.3 meV to the left of the upper mode- 
edge. Here phonon modification leads to reduction of SE 
rate, 7 < 7— (also Eig. 12(a)). Thus when phonons influ¬ 
ence the SE rate, the feeding of the waveguide mode-edge 
(solid line) is lower. 

The extent of feeding of the upper and the lower waveg- 
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FIG. 11: (Color online). Normalized projected spectra 
at T = 40 K in a CROW, for a QD (at 0 meV) located 1 meV 
to the right (a) and 0.3 meV to the left of the left of the upper 
modeedge (b). The orange (light) solid (black (dark) dashed) line 
plots spectra in case when phonons do ( do not ) influence the SE 
rate, 7 (7). The parameters are = 85 /reV, 70 = 7(i = 1-0 fieV. 


uide mode-edge by the QD can be quantified using peak 
intensity ratios as upper branch ratio, Ru = Uq/Iq and 
lower branch ratio, Rl = Lq/Iq- These intensity ra¬ 
tios is then plotted (Fig. 12(b)) as the QD is scanned 
in frequency across the waveguide LDOS at T = 40 K 
(see Fig. JT (b, c)). The blue (dark) and orange (light) 
solid lines plot the lower and upper branch ratios re¬ 
spectively, when phonons do not influence the SE rate 
( 7 ). The magenta (light) and black (dark) dashed lines 
plot the lower and upper branch ratios respectively, when 
phonons influence the SE rate, ( 7 ). The SE modification 
factor X = 7 / 7 — in presence of phonons as a function 
of wavguide LDOS frequency is also reproduced for ref¬ 
erence in Eig. 12(a). The broad side-bands in the peak 
ratios about the waveguide mode-edges appear due to 
intensity enhancement due to phonon side-bands. The 
sharp peak in the ratio happens near the mode-edge as 
the QD ZPL is scanned across the sharp waveguide mode- 
edge (width = 14 /ieV). The lower branch ratio Rl is 
stronger and wider than the upper branch ratio Ru inside 
the waveguide band, because lower (upper) branch ap¬ 
pears due to phonon emission (absorption). Eor the same 
reasons, the upper branch ratio Ru is larger and broader 
than the lower branch ratio Rl^ outside the waveguide 
band. In the case when phonons do not influence SE ( 7 ), 
as the QD is tuned across the wave-guide, stronger feed¬ 
ing is expected when the QD is aligned to the right of a 
mode-edge, due to stronger phonon emission—. This is 



FIG. 12: (Golor online), (a) Spontaneous emission enhancement 
factor X (at T = 40 K) plotted in loglO scale, as the frequency of the 
QD is scanned across the waveguide LDOS. (b) Lower and upper 
branch peak ratios plotted as the QD frequency is scanned across 
the waveguide LDOS at T = 40 K. The blue (dark) and orange 
(light) solid line plots lower and upper branch ratios respectively, 
when phonons do not influence SE rates (7). The magenta (light) 
and black (dark) dashed lines plot lower and upper branch ratios 
respectively, when phonons do influence SE rates (7). ouq marks 
the band-center frequency of the waveguide. The parameters are g 
= 85 /i-eV, 70 = 7 d = 1 A^eV. 


true for the lower branch ratio Ru without phonons (blue 
(dark) solid) which is stronger inside the waveguide band 
compared to outside. The reverse happens though for the 
upper branch ratio Ru without phonons (orange (light) 
solid). In case of the upper mode-edge, the reverse hap¬ 
pens, because the overall SE rate ( 7 , Eig. .9(a)) of the QD 
reduces as it moves outside the waveguide band. This 
effect is somehow suppressed when the phonon modifica¬ 
tion of SE rate (7 ^ 7 , Eig. 12(a)) is accounted for (black 
(dark) dashed). The feeding ratio is substantially larger 
in this case (almost twice), compared to the case where 
such effects are ignored ( 7 ). Phonons have an opposite ef¬ 
fect however at and near the mode-edge inside the waveg¬ 
uide band. Here the feeding is suppressed (dashed) due 
phonon induced reduction of SE rate^. 

A better understanding of the mode-edge feeding can 
be obtained by plotting the ratios as a function of temper¬ 
ature and detuning (Eig. 13(a)). Due to very small spec¬ 
tral overlap between the ratios (Eig. 12. (b)), a sum of the 
two ratios Ru Rl is plotted. The intensity of the side¬ 
bands increase with respect to the mode-edge peak with 
temperature, as the phonon sidebands become stronger. 
At very low temperatures, when phonon absorption is 
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negligiblei^, the feeding is stronger as expected, when the 
QD is aligned to the right of mode-edge. The trend how¬ 
ever reverses for the upper mode-edge as T is increased 
(12 (h), dark dashed lines). This is because the total SE 
rate inside the upper mode-edge is always greater than 
outside for the highest temperatures considered in the 
plot. 

The feeding ratios increase almost by a factor of 10 
(Fig. 13.(b)) when a temperature dependent pure dephas¬ 
ing term of the form 7 ^ = 1 -f0.95(T-l) /ieV— is added. 
Many of the finer features due to phonon modified spon¬ 
taneous emission 7 are however lost. The phonon modi¬ 
fied SE rate 7 is around 2 jaeV (T = 40 K) in the interest¬ 
ing regions outside the waveguide band, which is much 
smaller than 7 ^^. Thus the line-width of a QD is domi¬ 
nated by pure dephasing and do not change much as the 
QD frequency is scanned across the waveguide. Hence 
the intensity ratios are always stronger when the QD is 
to the right of a mode-edge, because phonon emission is 
stronger than absorption. 


IV. CONCLUSIONS 

In summary, we have provided an in-depth study of 
the incoherent emission spectra from QDs coupled to 
structured photonic reservoirs. We have carefully out¬ 
lined the approximations made in different theories and 
applied them under suitable conditions to study emission 
spectra from photonic crystal cavities and waveguides. A 
non-Markovian approach is developed using a correlation 
expansion treatment for coupled QD-cavity systems, and 
compared to the numerically much simpler approaches of 
cQED and bath polaron master equation theories, and 
we discuss the pros and cons of these approaches for 
modelling QD cavity systems. We also demonstrate the 
failure of the linear susceptibility approach for explain¬ 
ing far off-resonant cavity feeding. Using our general 
bath polaron master equation, a comprehensive study of 
the quantum emission spectra from weakly coupled slow- 
light waveguides is provided where effects of the non-local 
LDOS on the SE rate are clearly seen. We highlight that 
our general polaronic reservoir theory is completely gen¬ 
eral and can be applied to other structured reservoirs 
such as photonic band-edges, coupled cavities, finite-size 
waveguides and photonic molecules. 


Appendix A: Coupled-mode spectrum from an 
inverted atom or electron-hole pair 

In this appendix we compare the coupled mode spec¬ 
trum calculated using the polaron cQED ME (12.), for 
a QD excited by a weak incoherent pump 
solid light line. Fig. 14) and that obtained from initially 
excited QD (i.e (a+cr-) (t = 0) = 1) (5'™-"^™’^^^ dark 
dashed line. Fig. 14b The QD-cavity are off-resonant 
{Acx' = 1 meV and the spectras show perfect agree- 
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FIG. 13: (Color online). A sum of the feeding ratio Ru + Rl 
plotted as a function of temperature and detuning. The parameters 
are = 85 /i-eV, 70 = 1 jJ^eV and 7^^ = 1 jieV for (a) and 7^^ = 1 
+0.95(T-1) /xeV for (b). The color scale in (b) is saturated at a 
ratio of 8 to show features at lower temperatures. 


ment as long as the QD is weakly excited in the first 
case (5'™“^^^^). The reader should recall, that for an 
off-resonant high Q cavity (Fig. £.(b)), the coupled mode 
spectrum derived using the polaron cQED 

ME is close to the physically correct coupled mode spec¬ 
trum calculated using correlation expansion ap¬ 

proach. Thus the initial condition of an inverted atom 
correctly reproduces the spectra from a weakly driven 
QD. This clearly validates the initial condition used in 
linear susceptibility approach for calculation of linear 
spectra (Sec. IIP4b 

When a QD is initially excited, the system does not 
have any steady-state excitation (cavity/QD) in absence 
of a drive and the coupled mode spectra is calculated by 
modifing (421 ^ as. 
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For the calculations of Fig. lA (light solid line) the weak 
incoherent pump produces a steady state QD excitation 
of (cr+cr“) (t —> oo) ~ 0.01. When the pump is increased, 
the calculated spectra light solid line) starts 

differing from the inverted atom spectra 
dark dashed line) when (cr+cr“) (t ^ oo) > 0.1 (not 
shown here). This corresponds to a E/g > 0.01. 


FIG. 14: (Color online). Normalized emission spectra at T = 4 K 
for a coupled QD-cavity system under off-resonant condition. The 
thick light solid (dark-dashed) line denote coupled-mode spectra 

calculated using polaron cavity-QED 
ME approach for an incoherenty excited atom (initally inverted 
QD). The main parameters are g = 100 /xeV, k = 65 7 o = 5 

geV, = 55 ^aeV, = 1 meV. 
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